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Abstract 



The interaction between three localized vibrational modes is shown to be 
as relevant for the lifetimes of localized modes as the interaction involving 
two localized and one extended, and one localized and two extended modes. 
This contrasts with previous views. I support my arguments by a numerical 
study of a strongly disordered linear atomic chain. 
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Lack of a periodic structure makes the character of the vibrational states (VS) in glasses 
rather complex. Their spatio-temporal properties depend on the spectral region of interest. 
The lowest frequency VS are propagons (p, terminology from [|IJ) - sound waves with a 
wavevector (of magnitude 2ir/X), polarization and velocity. They elastically scatter off local 
structural imperfections, and as long as their mean free path £ > A, the scattering does 
not destroy their wave character. What happens at the Ioffe-Regel frequency ujjr, where 
£ ~ A (and the wavevector concept becomes invalid) is still a matter of debate. Since ujyr is 
very low (less than 5 meV for amorphous silicon (a-Si) where the maximum frequency is ~ 
80 meV), VS with uj > ujir are responsible for thermodynamic and transport properties of 
glasses at temperatures ranging from several Kelvins up to the melting point. 

Presently there are two different scenarios for what happens beyond the Ioffe-Regel limit, 
(i) The "fracton" model postulates that all VS with uj > c<j ir are "locons" (/), localized 
modes which are anderson-localized Q because of topological disorder. The relation to 
"fractal" geometry is an historical accident, (ii) In the "diffuson" model |||| the majority 
of modes are "diffusons" , extended and non-propagating modes, carrying energy diffusively. 
Numerical simulations on realistic models H0] show that locons form only a small portion 
(~ 3% in a-Si) of the highest-frequency modes. 

In harmonic approximation the VS will never equilibrate, but anharmonic interactions 
are always present. What is the rate of equilibration in glasses? Experiment [|J indicates 
that at a liquid He temperature in a-Si (a) high frequency {uj > 10 meV) VS live longer as 
frequency increases, and (b) the VS live on nanosecond timescales. In contrast, phonons in 
crystals decay faster as frequency increases (the number of decay channels increases), and 
lifetimes are picoseconds. Model (i) seems to agree with the experiment in the following 
way. Locons I can decay via three processes: (1) I — > p'p", (2) I — > p'l" or lp' — > I", and 
(3) I —>■ VI" or IV — > I". Process (1) is kinematically forbidden for high frequency locons 
with uji > 2co>ir. Processes (2) represent "propagon-assisted hopping" between locons I and 
I". Because uj p i {oji,uji"), uj\ ~ uji», and the locons / and I" spatially repel each other 
||. As the frequency and the inverse localization legth L" 1 increases, the overlap between 
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the locons decreases and so does the hopping rate. Processes (3) have been neglected by 
model (i) on the grounds of a small probability of three locon overlap. Hopping is therefore 
assumed to be the only locon decay process, and because it has property (a) and can be 
"adjusted" to fit (b), the experiment seems to be rationalized. 

Recent numerical calculations on models (i) and (ii) |1J have shown that decay properties 
of both diffusons and locons are similar to those of crystalline phonons. It has been argued 
U that the assumption of the negligibility of (3) is incorrect and is the reason behind 
the discrepancy between the numerical results and model (i) predictions. The experiment 
|| therefore remains to be explained. In this paper I show, using a one-dimensional (Id) 
numerical realization of model (i), that extended and localized modes decay on the same 
time scale and that the three locon interaction is dominant in the decay of locons with 
LU > 2u r . 



Consider a linear chain of N uniformly spaced atoms connected by random springs |]T0|.|TT 
Periodic boundary conditions are assumed. In terms of displacements u a («n+i = Ui) of 
atoms a from equilibrium, the potential energy of the system can be expressed as 

1 N 

V = -Y^K a {u a -u a+l ) 2 . (1) 

1 a=l 

Here K a =K (l +£ a ), where £ a are random numbers uniformly distributed in the interval 
<-b,b>, and K =10.6 eV/A 2 is chosen to simulate a hypothetical linear silicon chain (the 



silicon atom mass M) with maximum vibrational frequency uj ma , x = 2yK Q /M « 80 meV in 
the case of no disorder (b=0). The length scale is the interatomic spacing, the exact value of 
which is not important here. For the case N=3000, vibrational eigenstates and frequencies 
were found by exact numerical diagonalization. This number is sufficiently large to ensure 
that the conclusions will not be affected by finite size effects. The goal is to numerically 
realize the "fracton model" scenario with majority of VS localized. As shown later, the 
strong disorder value b=0.7 suits this purpose and will be used from now on. The "clean" 
case, b=0, will serve as a reference. 

Figure [l] shows vibrational density of states (DOS) for the above model and for the case 



b=0. The latter can be trivially solved by introducing wavevectors. Consequently its DOS 



has the 1/y^max — u 2 dependence and a Van Hove singularity at u max . As the disorder 
increases, the description of VS in terms of wavevectors is less and less valid. For b=0.7, 
the Van Hove singularity is washed out and DOS becomes nearly flat. 



In a Id infinite system any disorder causes all VS to be localized [fL^] , with a trivial 
exemption for the zero frequency mode corresponding to a translation of the system as a 
whole. However, if N is finite there are always low frequency modes with L larger than the 
system size. They appear to be extended, sound wave-like modes, although on intermediate 
timescales they behave as diffusons ||13|| . Only after the system size is increased, they become 
manifestly localized. A good measure of the number of atoms participating in the vibration 
of mode i is the participation ratio Pj. Its inverse is defined as 

N 

W = E(<) 4 , (2) 

a=l 

where e\ are normalized £ a (e^) 2 = 1) vibrational eigenstates with frequencies cjj. I plot 1/P 
in Fig. |2|. 1/P grows monotonously with frequency, the dependence being close to quadratic 
for u > 10 meV. At lower frequencies 1/P ~ 1/N, since L > system size (for infinite N it 
has been predicted jnj that 1/P ~ u 2 as a result of Id elastic sound wave scattering). To 



accurately locate the mobility edge u c I use the Thouless criterion |15| , that Au/8u should 
exhibit a sharp drop at u c . For a given mode, Au; is the locally averaged change of the mode 
frequency under the change of boundary conditions from periodic to antiperiodic, and Su 
is the local average level spacing. I find it convenient to consider the cumulative quantity 
t(uj) = Jq du' 'Au(u') I '8u(u') , which is plotted in the inset of Fig. [|. The mobility edge 
(vertical line) is u c « 9 meV and Au/8u at u c is less than 1%. Locons with u ^ u c form 
~ 90% of the spectrum, enough to simulate model (i). 

To get a qualitative understanding of how VS decay, consider a simplified version of 
zero-temperature decay rate formula (full formula see e.g. fT6|j ): 

fi = *52{Jijk) 2 8{ui - Uj - to k ), (3) 



where the sum is over all modes j and k, and 

Jijk = Va e a e a e a (4) 
a 

measures the overlap between modes i, j, and k. r] a is a binary random variable having 
values 1 or -1. Only "fission" decay processes (i — > jk) are allowed at zero-temperature. 
"Fusion" processes (ij — > k) appearing at finite temperatures will not be considered here, 
since they can be handled analogously. In Id VS have a partial memory of a wavevector, 
and therefore a partial phase coherence. The reason is that the n-th vibrational eigenstate 
has n-1 nodes (if the zero-frequency mode is n=l) jnj. rj a is inserted into Eq. |] to eliminate 
this coherence, so that the results will not be substantially affected by dimensionality (apart 
from the DOS effects). Imagine VS with amplitudes independent of frequency. Then an 
expression similar to Eq. |3] can be obtained for the decay of the VS caused by a set of 
nonlinear springs with potentials ~ i] n u' + O(w^), externally attached to each atom. 

The spectral dependence of T is shown in Fig. |. The case b=0, plotted in (b), does not 
represent the decay of Id crystalline phonons, since the phases of these are randomized by 
r) a in Eq. f|. It, however, shows the hypothetical behavior of random-phase extended modes, 
such as occur below the mobility edge in 3d. The cu-variation of T follows the joint density 
of states (JDOS), J2jk ~ u j ~^k), and provides a benchmark when discussing the decay 
behavior of locons. Data in (a) at u u c are much more scattered than in (b), since the 
locon decay involves some overlap statistics (explained below). Therefore local averages are 
used to represent the data. V in (a) increases with frequency and follows JDOS, having 
similar magnitudes as in (b). According to model (i) there would be a sharp decrease in V 
beyond 2uo c . This certainly does not happen. 

The reason why the behavior of V differs from the one predicted by model (i) is 
that processes (3) can not be neglected. In fact they are dominant decay processes 
for high frequency locons. Consider a decay i — > jk, where i, j, and k are locons 
with corresponding participation numbers Pi, Pj, and Eq. [3] can be rewritten as 

Ti = J J dujdiJkN(uj)N(uJk)Ji(uj,Uk)5(u)i — ujj — u>k), where N(ujj) is the total number 



of states at u>j, and Jf(u>j,ujk) is Jfj k averaged over frequency shells of Uj and u>k- A simple 
scaling argument ]I| shows that Jf(ujj,uj k ) does not explicitly depend on P^ Pj, or P k . It 
goes as follows. The maximum decay magnitude among Jfj k is found when i lies entirely 
within an overlap of j and k, and scales as 1/PjPk (e* ~ 1/y/Pi, etc., and the sum in Eq. |] 
is over random-sign numbers and scales as \fP~i). The probability for this overlap to occur 
scales as PjP k /N 2 . An important assumption used here is that locons i, j, and k, having 
different frequencies, are spatially uncorrelated (illustrated below). The average Jf(uj,u k ) 
therefore scales as 1/N 2 and Tj is N-independent (and follows the JDOS). The small prob- 
ability of the three-locon overlap is exactly compensated by the large magnitudes of decay 
matrix elements, when the overlap occurs. As for hopping processes I — > p'l", one has to 
distinguish between two cases. If <C uji-,uJi» (as in model (i) and this paper), the two 
locons / and I" spatially repel each other || and the decay rates decrease as L of I and I" 
decreases. On the other hand, if u p > is not small, one can apply the above scaling arguments, 
since majority of I and I" are not spatially correlated. Decay rates of this "uncorrelated" 
hopping are therefore independent of the localization character of VS involved and follow 
the JDOS. When a low frequency (u>i < 2u c ) locon decays into two propagons, I — > p'p", 
the overlap probability is one (Pj,P k = N) and the case is trivial, similar in behavior to 
three locon and "uncorrelated" hopping processes. These arguments do not depend on d 
and worked, as far as one could tell, in 3d also 

Figure | illustrates the foregoing discussion. Data are represented by running averages. 
Decay I — > VI" in (a), obtained by considering only j and k with ujj, uo k > uo c in Eq. |3], 
dominates the high frequency region and grows monotonously with frequency. When, on 
the other hand, restricting, say modes j to be propagons (u>j < u c ) , one gets hopping I p'l" 
which decreases as frequency is increased. Since u c is large enough to allow hopping between 
locons from different spectral regions, the decrease is quite moderate. The graph (b) in Fig. 
f| shows hopping when propagon frequency is restricted by u p < 1 meV. Computed rates 
are scattered over several magnitudes. In a sample of decay rates of locons confined into a 
small frequency band (of the width, say, 1 meV), majority have negligible values, but there 
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are ^ 10% of large magnitude rates giving rise to large local averages (they come from a 
residual overlap between locons with similar frequencies). As the graph shows, local medians 
sharply decrease with frequency and deviate strongly from local averages. Since the local 
distribution of decay rates is found to have a strong central tendency around small values 
(the area in the tails is small), the medians better represent computed data. For other decay 
types local medians and averages are similar. 

Finally, I present figure || to support the argument that locons from different spectral 
regions are spatially uncorrelated, as opposed to locons with similar frequencies, which repel 
each other. N=1000 is used for this purpose. The spectrum is divided into 5=1 meV wide 
bands, each band containing 10 VS on the average. For every mode % in a band only atoms 
a having (e % a ) 2 > 0.2 (e^^) 2 are shown, e % max being the largest magnitude found among the 
coordinates of the vector e\ Localization starts somewhere beyond 15 meV, where fewer 
than P Pd 200 atoms participate in the vibration. There is clearly no statistically relevant 
overlap between the locons lying within a frequency band. On the other hand, locons from 
different spectral regions show no sign of repulsion, their spatial positions are essentially 
uncorrelated. 

In summary, I have used a Id random spring atomic model to demonstrate that the inter- 
action including three locons is important for the locon decay processes. I have shown that 
extended mode and locon decay rates have similar magnitudes, increasing with increasing 
frequency and closely following JDOS. Hopping decay rates decrease as frequency increases, 
the decrease is sharper as the frequency difference of hopping locons gets smaller. 

I thank P. B. Allen for useful comments on the manuscript, and J. L. Feldman and S. 
Bickham for stimulating discussions. This work was supported by NSF Grant No. DMR 
9417755. 
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FIG. 1. Calculated DOS for the 3000 atom linear chain with the spring constant disorder of 
b=0.7 (solid line) and b=0 (dashed line). 



FIG. 2. Inverse participation ratio 1/P and the cumulative Thouless number r (inset) as a 
function of frequency. The vertical line in the inset is the mobility edge uj c ~ 9 meV. 



FIG. 3. Calculated T (circles represent uniformly sampled 1/10-th of data) for b=0.7 (a) and 
b=0 (b). Scattered data in (a) are represented by their running averages (solid line). Vertical lines 
indicate uj c and 2u c , and dashed lines are JDOS. Units are arbitrary, though the same in (a) and 
(b). 



FIG. 4. Calculated f (a) for the three locon (solid) and hopping (dashed) decay processes, and 
(b) for hopping between locons from the same frequency region (circles - see caption to Fig. |||). 
Vertical line in (a) indicates io c . Units are the same as in Fig. |3[ 



FIG. 5. Localization of VS in the 1000 atom linear random spring model. For each mode only 
the atoms vibrating with amplitudes higher than a certain value (see text) are shown. The solid 
line is the particiation ratio P. 
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